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An single Israel layer can be created when two metrics adjoin with no continuous 
metric derivative across the boundary. The properties of the layer depend only on 
the two metrics it separates. By using a fractional derivative match, a family of 
Israel layers can be created between the same two metrics. The family is indexed by 
the order of the fractional derivative. The method is applied to Tolman IV and V 
interiors and a Schwarzschild vacuum exterior. The method creates new ranges of 
modeling parameters for fluid spheres. A thin shell analysis clarifies pressure/tension 
in the family of boundary layers. 
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I. INTRODUCTION 



There is long standing interest in fluid sphere solutions, largely because of their astrophys- 
ical implications. An astrophysical model is often an interior fluid sphere metric matched 
to a Schwarzschild vacuum or Kottler exterior across a bounding surface. The standard 
technique matches metric functions and extrinsic curvatures on the boundary. When the 
extrinsic curvatures do not match, an Israel boundary layer can be created. The layer 

depends only the properties of the two bounding metrics. Methods that will create a fam- 
ily of surface layers between the bounds could prove useful in exploring models of spheres 
with variable crusts. One way of creating variable surface layers is to modify the boundary 
conditions at the fluid- vacuum interface. 

While an extrinsic curvature match is the boundary condition currently most used, there 
are three types of boundary conditions that have been used to match analytic solutions 
across non-null boundaries. The three methods have been discussed by Bonnor and Vickers 
a. and they al, involve deHva.ve= of .he n,e.nc funct... The bounda.. condMo. can be 
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generalized by broadening the idea of derivatives to include fractional derivatives 

There are two simple ways to proceed with the generalization. The first is to assume a 
straight fractional derivative match on the boundary metrics and then to use the fractional 
relations in the usual formalism for the boundary stress-energy. This would be a generaliza- 
tion of the Lichnerowicz boundary condition. It would not generalize the extrinsic curvature 
to fractional values. The second would be to use fractional derivatives to define a fractional 
extrinsic curvature and then use it to define a fractional boundary layer. This would be a 
generalization of the usual Lie derivative to fractional values. The use of fractional calculus 
is motivated by the possible fractional nature of the growth processes forming the boundary 
layer. Fractional transport processes are one of the main areas of application for fractional 
calculus, and boundary layers formed by these processes could reflect this fractional forma- 
tion process. 

Beyond the fractional generalizations of techniques and tensor functions, one must con- 
sider the various definitions of fractional differentiation. Use of fractional calculus in diverse 
areas of physics has increased enormously since fractional derivatives were first considered by 
Leibnitz and L'Hospital in 1695. Many different definitions have been proposed for differ- 
ent applications. In this article we use the Caputo form of the Riemann-Liouville and Weyl 
definitions. The Caputo derivative is an integral transform of the regular partial derivative 
and preserves zero fractional derivative of a constant. While considering generalizations of 
relativistic gravity to include fractional calculus, the different definitions must be explored 
to determine their applicability. 

This work has two goals: first to develop a variable layer model that could be applied to 
astrophysical problems, and second to better understand the role that fractional derivatives 
might play within a general relativistic framework. In this article we will apply the first 
method and use fractional derivatives to create a family of Israel boundary layers between 
two bounding metrics. The family is parameterized by the order of the fractional derivative 
and may be used to model fluid spheres with variable crusts. Even when a regular derivative 
match is possible, the fractional match will broaden the parameter ranges for the fluid 
interior. 

In the next section we discuss the metrics and describe the boundary layer. In the 
third section several models are considered: the Misner-Zapolsky (MZ) solution , and 
Tolman's solutions IV and V Q]. The thin shell pressure balance is treated in the fourth 
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section and some details of the fractional match are discussed in section five. Details of the 
fractional derivatives and the standard fluid sphere formalism are given in Appendices. 

II. THEORETICAL FRAMEWORK 
A. The spacetime 

The two regions to be considered are covered by an exterior Schwarzschild solution bound- 
ing an interior spherical fluid. The metrics are, with functions ipsch = 1 ~ 2mo/?/, i^{r), A(r), 
H{r) 

Exterior: glfdx'^dx'^ = -ipschdt^ + i^sch^y^ + V^^n^ (1) 
Interior: glf'^dx^'dx'^ = -e^dr^ + e^dr^ + H^dQ"^ (2) 

The bounding surface is located aX y = in the exterior and r = Rq in the interior. The 
corresponding normals to the surface are 

Exterior: n^dx^ = 'ipsch'^^y (^) 
Interior: ri^^dx'^ = e^^'^dr (4) 

Fractional derivatives leave these metrics unchanged. Our fractional extension provides a 
crust layer between the interior and exterior metrics. 

B. Matching Conditions 

On the boundary, the metric match conditions are 

(1 - 2mo/i?o) = e^(^«) (5) 
Ro = H{Ro). 

The second matching condition is the extrinsic curvature match, K^, on the bounding sur- 
face. If the curvatures do not match, an Israel boundary layer is created. The stress-energy 
content of the Israel layer is constructed from the mismatch in the extrinsic curvatures. The 
stress-energy of the boundary layer is P| 

-87rS^= <K-,>-<K>gl (6) 
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Here K = K^. The stress-energy components on the boundary are: 



-87r5° = [<K',>-< + 2K'g > g',] = -2 < K', > 

-SnS', = -SttS; = [< K', >-<K + 2^.' > 9e] = "[< Ke > + < Ko >] 



and the stress-energy of the boundary is 

1 
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C. Match of fractional derivatives 



The stress-energy of the Israel layer is evaluated on the boundary between the interior 
and exterior metrics. The actual finite thickness boundary layer is modeled by the single 
bounding surface at r = i?o- The stress-energy content is governed by regular derivatives of 
the metric functions. The metric match coupled with some derivative match of the metric on 
the layer, sets relations between the parameters of the interior and exterior solutions. With 
the usual extrinsic curvature or other derivative matches, the properties of the layer are set by 
the parameters of the bounding metric. With a fractional match, the order of the fractional 
derivative enters along with the other parameters and a family of fractional boundary layers 
is created. The fluid sphere examples considered in this paper have boundary metrics of the 
form 

ds^ = -F{r)de + r^d^l^. 

The fractional match is applied only to the differing part of the Israel layer metric, the qqq 
metric potential. The actual calculation of the fractional derivatives involves a choice of 
definition. We use the Caputo definition (see Appendix A) with the (0 < r < Rq) Riemann- 
Liouville limits for the interior and the [Rq < r < oo) Weyl limits for the exterior. The 
limits themselves, as well as the choice of different limits for interior and exterior derivatives 
reflect the non-locality of the fractional derivative operation. Non-locality in fractional time 
derivatives is an expression of system memory |9j] . It has proven especially useful in modeling 
jump processes with long wait times Q]. Similarly, spatial non- locality implies the derivative 
on the boundary depends on values away from the boundary; fractional spatial derivatives 



have been useful in modeling processes with very large jump distances When the jump 
distance depends on the jump time, fractional time and spatial derivatives enter into the 

are static but the structure of the 
boundary layer could reflect the transport process. The fractional matching condition is 

-Ro oo 

1 f(P^FAx)ldx'^ , (-1)""^ rd'^FE{x)ldx'^ , 
dx = — : / — ^ dx 



T{n-a)J (r-a;)"-"+i T{n-a 

Ro 

and is applied at r = Rq. We note that the single layer at r = Rq only approximates a 
boundary of finite thickness and that using a non-local operator might be a better approx- 
imation to the actual match over a finite thickness than the usual derivative match over a 
zero thickness surface. 

In the next sections, we apply the formalism to Tolman IV and V solutions. 

III. MODEL CALCULATIONS 
A. Tolman's Solution V 

1. The solution 

We consider a parametrization of Tolman's Vth solution The metric, with con- 

stants n and C, is 

ds^ = -{r/rof'dt^ + a(l - aCr^+')~^dr^ + r^dVl^. (8) 
The parameters formed from n are 

Ni = 4:n/{l + n), A^2 = 1 + 6ra + 

N2 N,{l-n) 
(l + n)2' (l + 3n) ■ 

The interior density and pressure for this solution are 

An 1 

8'rp=(^)^ + C(3 + (,),•', (9) 

^ N2 1 + n 

For C = 0, the solution reduces to the MZ solution This solution was 

originally used to describe neutron star models with equation of state P = np. The solution 
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with C = does not admit a zero-pressure boundary, the C 7^ solution does. Both solutions 
are singular at the origin and are used to represent an ultrahigh density core. The MZ 
solution, lacking a vacuum boundary, is generally matched to a gaseous envelope. It may 
be more realistic in some cases to match these solutions to a crust with surface stresses. 
For C 7^ 0, the zero-pressure boundary, i?^, relates constants C and n 

Ain?{l + n) 1 
[l + bn)N2'^^ 

Substituting for C, the pressure and density can be written as 



c = ^f^^i^^. (10) 



SirPc^o = ^ I -1 - -I- 1 (11) 





(r2 . 






4n / 




iV2 \ 



— ' / - > (n + 3) / , 



The condition Pc^o > requires 



Ro < Rz (13) 



Fractional matching will allow a broader family of sphere sizes. Below, we graph values 
for the case n = 1/3. For C = 0, there is no zero-pressure boundary and no constraint. 



2. Matching Conditions 

The matching conditions are the same for any C value. Matching the interior metric to 
vacuum Schwarzschild we find 



1 - 2mo/i?o = (i?o/ro)''^ 



[recall iVi = 4n/(l + n)] 



Rl+" ' ' ro' 'R^'T{Ni + l-a)' 
Combining the two relations we find the scaled radius of the interior is 

Ro ^ ^ r(i + a)r(iVi + i-«) ^^^^ 



2mo r(l + Ni] 

Note that the boundary radius is always greater than 2mQ. For a < 1, there are no limits 
imposed by Eq. ()14|) . For a > 1 we require 



1 + A^i > a 
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The metric parameter tq is described by 

= 1 + 



r(i + iVi) 



r(l + Q;)r(A^i + 1 - a) 
The sizes of the fractional spheres are discussed in section V. 



(15) 



3. The Crust Stress-Energy 



The stress-energy of the crust for general C is, with 70 := 1^1 — 2mo/ Rq 



Kq 



SttS^^ = SnSt 



70 - (1 + n)N;'/'^l-aCRl+' 

3 + 1/70 - (1 + n)(2 + N^)N-^'^^l-aCRl-'' 



2i?n 



For C 7^ 0, the boundary layer has a stress-energy content, [recall A^i = 4n/(l -|- n), N2 
l-\-'on-\- V?] 



8nS"^ = 8nSl = 



7o-(l + n)iV2 
1 



1/2 



1 - 



2Ra 



(1 + 5n) 
7o + l/7o-2(l + 3n)iV2-'/' 



(i?o/i?.)2+'' 



(16) 



nNi 



(1 + 5n) 



(i?o/i?.)2+'' 



For C = the fluid energy density and stress are 

871^0° = (2/^0) [70 - (1 + n)N;'^'' 

= (2/i?o) f(i?o/ro)'"/('+") - (1 + n)iV, 



(17) 



-1/2 



Stt^^^ = Stt^J = (l/i?o) [(1 - mo/Ro)ho - (1 + 3n)7V-i/2 
and describe a much richer modeling environment. 



B. Tolman's Solution IV 

1. Metric and Stress-Energy 

This solution describes an object with finite central presssure and density. A stiff fluid 
core is not possible in this model. The interior metric for this solution is, with constants A, 
B, and C 

l + 2rVA2 



-dr' + r^dn^ 



(l-r7C2)(l + 7-2/^2) 



(18) 



The interior density and pressure are 



l + 3(AVC2 + rVC2) l-rVC2 



l + 2rVA2 (l + 2rV/l2)2 

1 _ + 3r7C2) 



l + 2rV>12 

Constants A and C can be expressed in terms of the central fluid values. We have 



87t{pc/3 + Pc) 
2 



87r(pc/3 - Pc) 

Note that the central fluid equation of state (EOS) is constrained: Pc < Pc/^- The 
pressure boundary that occurs in the regular derivative match has size 



2. Metric Match 



The match to vacuum Schwarzschild provides 



B'{1 + nil A') = 1 - 2mo/i?o- 



The fractional match is 



^2 



mo 



>12 i?^-2 r(3 - a) 

B'^Rl = A^mo r(3 - a)r(l + a) 



Combining with the metric match we obtain 



^2 ^ ^2 (^o/mo) - [2 + r(3 - a)r(l + a) 



B 



2 ^ ^ 
Rq 



mo 



r(3-a)r(l + a) 
« _[2 + r(3-a)r(l + a)] 



i?o/mo > 2 + r(3 - a)T{l + a) 
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The boundary size depends on the central EOS as well as the order of the fractional 
derivative. 



mo 



mo 



r(3-a)r(l + a) 
47rm2(P, + p,/3) 







(25) 



3. Crust Stress-Energy 



We introduce scaled parameters r\ := RyA^, := RyC^, and := RyA^. 



8nS^ 



1 



1 



2 
2 

Rq 

1 

Rq 

+ 



1 - 2mo/i?o - 



\ l-rl){l + r\) 
l + 2r\ 



l + r\ 



B 



A ll + 3rl-r\ 



C 



l + 2r\ 



«)^/^ 9fe ialrY^' die 



+ 



(26) 



1 fi'lX),?/ 1 fi'oo,; 



l + r\ 



^_A ll + 2>rl-r\ 



l + 2r\ 



(27) 



B 



l + 2r2 



Some examples of radius variation and crust stress energy are given in Section V. 



IV. EQUILIBRIUM IN THE PRESENCE OF SURFACE STRESSES 

A. Stress-Energy 

The Israel layer is the zero-thickness idealization of a bounding layer with finite thick- 
ness, d. The physical crust runs from an outer boundary Rq to an interior fluid boundary 
Ri with d = Rq — Ri. We know the interior fluid solutions will satisfy the 0] Tolman- 
Oppenheimer-Volkov (TOV) equation. The Israel layers generated in this work are obtained 
by introducing a discontinuity in the derivative of Qoq. The analog of the TOV equation for 
the layer, requiring that the solutions remain static, will provide relations among the model 
parameters. To develop the TOV analog for the layer, consider the general static spherical 
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metric for an interior fluid with pressure P and density p 

ds"^ = -e'^dt^ + e^dr^ + r^dQ\ (28) 

The details of the field equations are given in Appendix B. The covariant derivative of the 
general energy-momentum tensor provides the conservation equation 

^ i^ + -)Tr+OT^ + eT^ = (29) 



dr '2 ' '2' W 

For an isotropic fluid matched to vacuum, this is the usual TOV equation 

dP u' 

^+2(^ + P)=0- (30) 
It is the analog of this equation which we want for the Israel layer. 

B. The conservation equation over a limiting shell 

Consider a bounding shell which will approximate a thin surface layer. The central radius 
of the shell is R with the outer boundary i?*-^-* = R + d/2, and the inner interior fluid 
boundary at = R — d/2. d is the coordinate shell thickness. In the c? — >■ limit, 

R ^ Ro- A general stress-energy Tj can be related to a surface stress-energy by 

T; = 5{l)Ste\e] (31) 

where ej^ is a tangent vector to the shell, and / the proper distance along a radial geodesic, 
I — Q^/'^d. The shell stress-energy has a perfect fluid analogue 

S'^ = aWW + T{h'^ + WW) 
n' = (0,e-^,0,0) 

where Sq/c^ = — (T(gm/cm^) and Sg = ±r(dynes/cm). Following Poisson we take / = 
on the hypersurface defined by R, with I negative for r < R and positive on the vacuum 
side, r > R. The TJ^ content of the shell can be described using a Heaviside function, 0(/), 
as 

t; = e{i)T,^+> + e{-i)Tt> + 6{i)s;. (32) 



11 



The last term will be zero for the 2 + 1 shell stress energy. Forming the derivative needed in 
the conservation equation we have 

or dr or dr or 

In the / — * limit we have 

frpr J] 

= - lim[5(/)— T,(-)n = - lim[(5(/)Pe^/2] (33) 
dr 1^0^ dr 1^0 

where the first term is zero with no radial pressure on the outer boundary. The stress- 
energy function evaluated at the inner boundary is P and A*-"-* is an interior metric function. 
Substituting into the conservation equation in the / —>■ hypersurface limit we have 

Pe^'-'/2 + (^Y2) 5° + (2/i?o) 5',^ = 0. (34) 



C. Evaluating (?rgoo 

The derivative, z/', on the hypersurface can be written as a difference equation 

_ iy{R + d/2)-iy{R-d/2) _ ujR + d/2) - ujR) + u{R) - u{R - d/2) 
' ~ {R + d/2) -{R- d/2) {R + d/2) -{R- d/2) 

_ p{R + d/2) - u{R) iy{R) - ujR - d/2) 
~ d ^ d 

Expanding, we can write 

iy{R ± d/2) = iy{R) ± i^'iR^^^]"^ + ... 
Substituting in the thin shell limit we have 

The first term follows from the Schwarzschild metric match and the second is given in 
Appendix B. We have 

u'{R) ^ {2m /R^ + 47ri?P)(l - 2m/ R)^^ 

where we have identified the Schwarzschild mass parameter with the interior mass of the 
fluid. Substituting into Eq.(j2D) we have 

-P(l - 2m/ R)-^/^ = {m/R^ + 2nRP){l - 2m/R)-^S^ + {2/R)Sl (35) 
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which describes the thin shell pressure balance. The classical limit of this equation follows 
from c ^ oo and is 

m , 2 , , 

If the fluid pressure at the interior boundary dominates, this can be interpreted as a tension 
in the shell balancing the outward interior fluid pressure at the boundary minus the inward 
pressure due to the gravitational attraction of the shell by the interior fluid. If the shell 
mass term dominates, the stress in the boundary layer will be a pressure. In the next section 
we explore the stress-energy structure of the boundary layer and will see parameter ranges 
with both layer tension and pressure. 



V. DETAILS OF THE FRACTIONAL MATCH 



A. Sphere Radii 

The sizes of the sphere are described by 

Tolman V : ^ = 1 + £ii±4£<Mi±^^li^ (36) 

2mo r(l + 4n/(l + n)) ^ ^ 

The scaled boundary radius, Ro/rriQ, for Tolman V is plotted as a function of alpha for 
various n in Figure 1. The overall effect of the fractional match is to increase the range of 
sphere sizes for a given EOS. The largest differences are for low n fluids where a very much 
smaller sphere radius is possible than for the zero-pressure match. The Buchdahl bound 
limits the ratio of 2mo/-Ro for fluid spheres whose Qqq component is continuous across the 
boundary and whose density is decreasing outward. We have matched fractional derivatives 
rather than first derivatives and it is not clear that the conditions of the Buchdahl bound 
are satisfied, but from Figure 1, it is seen that the Buchdahl bound, 21710/ Rq ^ 8/9; is not 
violated. 

The radius for Tolman IV is a cubic root of Eq. (j37p but some general description can be 
given. The modeling term in the equation is the denominator of the last term. Consider 
the factor 

ci ~ Anml pc/3 
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FIG. 1: Scaled radius vs fractional order 

describing an object with mass rriQ = Ntuq = N{2 x lO'^^A;^'), and central density and 
pressure of neutron star order, pc ~ 10^"^ kg /m^, Pc ~ 10^^ Newton/m^. Numerical scaled 
radius values using these values are described in Table I for a range of N and alpha values. 
For ~ 100 or larger, the last term in the cubic is negligible and the radius is essentially 
given by the limiting value 

Ro/mo ~ 2 + r(3 - a)r(l + a) 

The masses for these radii are well out of the neutron star range. The reflection symmetry 
about a = 1 is the result of a product equivalence of the two gamma functions for paired 
alpha values, i.e. a = (0.6, 1.4) give the same gamma function product. From Table I it 
is clear that the low values have masses and radii of neutron star orders of magnitude 
jlG^ . For example, for = 1, the radii are approximately Rq ~ 11.44mo ~ 17km . Smaller 
central densities, describing more ordinary fluid objects, result in much larger fluid spheres. 
For a central density of pc ~ lO^^kg/m^, the radii for a = 1, Ro{N) are -Ro(l) = 2228. 33mo, 
i?o(10) = 480.87mo, i?o(100) = 104.39mo, /2o(1000) = 23.32mo, with the values for larger 
and smaller a paired and increasing, just as for the larger central density. 
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N = 1 


iV = 10 


N = 100 


= 1000 


0.2 


13.24 


4.41 


3.55 


3.54 


0.4 


12.39 


4.10 


3.28 


3.27 


0.6 


11.85 


3.91 


3.12 


3.11 


0.8 


11.54 


3.81 


3.04 


3.03 


1 


11.44 


3.78 


3.01 


3 


1.2 


11.54 


3.81 


3.04 


3.03 


1.4 


11.85 


3.91 


3.12 


3.11 


1.6 


12.39 


4.10. 


3.28 


3.27 



Table I: Rq/Mo - Tolman IV-p^ ~ 10^^ kg /m^ 

B. Crust Stress-Energy 




FIG. 2: Ds = — STrmo^o vs fractional order 



Figures 2 and 3 describe the variation of the boundary energy density and pressure, in 
Tolman V, n = 1/3, as the size of the fluid sphere varies. Over a large part of the Rq/Rz 
range, a crust tension contains the interior fluid, with the tension increasing in size as the 
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FIG. 3: Ps = SvrmoS^ vs fractional order 

fluid sphere becomes smaUer, essentially acting to squeeze the fluid into smaller volumes. 
The results are similar for C = 0. 

For Tolman IV, the crust energy density is 



Sttct 



2 

Rq 



The modeling factor of importance is the term 



'i l-Rl/C^){l + Rl/A^) 
1 + 2Rl/A^ 



r-——— [1 - 47ri?gpe(l/3 - n)][(l + 47rp,i?g(l/3 + n) 

-VI - ^mo/Ro + \ 2q n /Q I — N 

V 1 + i?n87rpc(l/3 + n) 



A -'•'0 2 Q -'^0 



For nuclear central densities, pc ~ 10 g/ cm , this is 

20.35A^' X 10-^^. 



m, 
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In order to have real values we require 

N =1, i?o - 12mo, [0.293 x (1/3 - n)] < 1 
N = 10, Ro ~ 4mo, [3.26 x (1/3 - n)] < 1 
N = 100, Ro ~ 3mo, [183 x (1/3 - n)] < 1 

It is clear that the broadest range of central equations of state for nuclear central densities 
is for the lower mass objects. Higher mass objects require a central EOS very close to the 
1/3 limit. For smaller values of the central density the central EOS range is much broader. 
For a central density of pc ~ lOi^g/cm^, the modeling factor is 

mt m^o 

and for real values require 

:228mo, [0.10 x 10-^(1/3 - n)] < 1 

(on™ rn A'7 in— 2/1 /o „m ^ i 





1, Ro ' 




10, Ro 


N = 


100, Ro 




1000, Ro 



VI. CONCLUSION 



In this work we have examined a family of boundary layers created by matching frac- 
tional derivatives across a boundary. The boundary layers considered have structure which 
depends on the order of the fractional derivative. One of the reasons that fractional calcu- 
lus may be important for boundary layers is the mechanism by which a boundary layer is 
formed. One of the possible ways to build a variable density crust is by a diffusive process; a 
process whose underlying cause is Brownian motion. This motion, as analyzed by statisti- 
cal mechanics, involves diffusion, dissipation, and the fluctuation-dissipation theorem. The 
dynamical model of Brownian motion was provided by Langcvin in 1908 using a stochastic 
differential equation. It seems apparent from the nature of randomness that such macro- 
scopic stochastic equations are incompatible with the continuous and differentiable character 
of microscopic Hamiltonian dynamics. (Think of the conventional diffusion equation, with 
the diffusion process described by a second order spatial derivative.) Therefore, the math- 
ematical description rests on either ordinary analytical functions describing the dynamics. 
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or on conventional differential operators describing the phase space evolution. The differen- 
tiable nature of the macroscopic picture is, in a sense, a natural consequence of microscopic 
randomness. This means that use can be made of ordinary differential calculations on the 
macroscopic scale, even if the microscopic dynamics are incompatible with ordinary calculus 
methods. On the other hand, in the case where a timescale separation between macroscopic 
and microscopic levels of description does not exist, the non-differentiable nature of the 
microscopic dynamics is transmitted to the macroscopic level. Since fractional calculus has 



been shown to provide a good description for a range of diffusive processes 19|, one might 



expect a boundary condition based on fractional calculus would reflect the fractional growth 
process. An example, given by Allegrini, Grigolini, and West [2^, shows that a diffusion 
process generated by a fluctuation with no time scale at the macroscopic level generates a 
diffusion process well described by a fractional Laplacean. 

While diffusion is a possible mechanism for generating a layer with structure, the method 
of generating the fractional family of layers is independent of the production mechanism 
and generates a family whose stress-energy and size depends on the order of the fractional 
derivative. The result is a much broader range in fluid sphere properties. For example, 
in Tolman V, a much larger range of spheres sizes can be described with the fractional 
layer than without, with the energy density of the layer decreasing as the size of the sphere 
increases. The C ^ Tolman V spheres have a zero-pressure boundary solution. For spheres 
smaller than the zero-pressure sphere, the layer has a tension, while for spheres larger than 
the zero-pressure sphere, the layer has positive stress over much of the range of the fractional 
order. The fractional boundary could prove to be a valuable modeling tool in more realistic 
neutron star models. 

The range of stress-energys in the fractional boundary layers implies differences in struc- 
ture as a function of the fractional order. The differences in density could be modeled in 
several ways: for example, with different crust materials or different incomplete fluid cover- 
ings (tilings) [2^. The layer itself is a model of a thin crust and there could be differences 
in the geometry of the 2+1 shells that fill the real crust. The interior fluid geometry and the 
exterior Schwarzschild vacuum do not have to match for a crust with finite thickness. 

The models presented in this paper matched an integral transform of the regular derivative 
across a spatial boundary. It is not a fractional generalization of general relativity, but a 
fractional genralization of a boundary condition. The next step is to explore the range of 
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fractional generalizations of other sets of boundary conditions and other derivative definitions 
to applicable spacetimes. 

APPENDIX A: DERIVATIVES 
1. Regular 

For functions / and F continuous on [a,b] e Reals 

F(x) - Jf{t)dt 

a 

F{x) is differentiable such that dF/dx = f. The n*'* integer derivative is simply d^F/dx^. 
For example 



dx" {k-n)\ 

With the gamma function this is 



dx"" T{a-n + 1) 
The gamma function, r{z), is defined as 

oo 

T{z) = Je-H'-^dt, r(i/2) = ^AF, r(i) = 1 



r(n + 1) = nr{n) = n\ for n > 0. 

2. Fractional 

a. Riemann-Liouville 

The Riemann-Liouville definition for the a fractional derivative of f{x) is, with a > 0, 

d" 1 (i" f(f] 

D"f{x) = -^fix) := A^^cit (Al) 

■'^ ^ dx'^-'^ ' T{n-a)dx'^i{x-tY-''+^ ^ ' 

where n is the smallest integer larger than a when it is fractional, that is n = [a] + 1. In 
the a = 1 limit, the derivative produces the integer result. The constant c in the limit of 
the integral is usually set to (Riemann definition) or to — oo (Liouville definition). For 
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example, the Riemann-Liouville derivative of x'^ for a < 1 with n = 1 we have 

' 

-rt'=a;-"(l- -)-" dt 

■ 

Using the definition of the beta function 



1 


d 


r(i- 


(x) dx 


1 


d 


r(i- 


a) dx 


1 


d 


r(i- 


a) dx 



r(p + q) 



we have 



^ 1 dx-"+i+'^r(A; + l)r(l-a) 



D'^x' 



r(l - a) dx Y{k + 2 - a) 
(1 + A; - ol)x^ 



.-a nk + i) 



T{k + 2-a) 

_^k-a r(fc + i) 
T{k + l-a) 

For q; = 1, this is the usual result D^x'' — dx^~^. Note that, for k = —1, this operation fails. 
This is an example of one of the problems encountered in applying fractional derivatives 
to general relativity. Not all definitions of fractional derivatives work for all functions. The 
fractional derivative of x^ for a > 1 is identical to the fractional derivative for a < 1. For 
this function, the derivative is continuous across the a — 1 boundary. 

D"^^ = J- jljtV-"(i - -f-'^ dt 



r(2 — a) dx'^Q x' 



r(2 — a) dx"^ 

k-a r(i + fc) 



X 



r{l + k-a) 

One should note that the Riemann-Liouville fractional derivative of a constant is not zero. 



b. Caputo 



The Caputo derivative is the integral transform of the regular derivative and is found by 
moving the derivative in the Riemann-Liouville definition inside the integral to act on the 
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function. We have 



D-f(x) = — f , '''",', dt (A2) 



Example: f{x)^x\ 13>Q 



1 X d^fk 

D"{x - a)' = J — -dt 



For a < 1, n = 1 we have 



1 X 

D«(a; - af = — Jkt^'\x - ty^dt 

-jw''-\l-w)-''dw 



r(i - a) 







which is identical to the Riemann-Liouville derivative for this function. This derivative also 
is not defined for for A; = —1. In general relativity, one of the spacetimes one would like 
to treat is vacuum Schwarzschild but the Riemann-Liouville derivative will not give finite 
answers for the 1/r structure. The derivative of 1/r can be taken with the Weyl derivative. 



c. Weyl 

The Weyl derivative differs from the Riemann-Liouville derivatives over the range of the 
fractional transform. To take the fractional derivatives of 1/r we use the Weyl derivative 
over the range {Rq,oo). The Weyl derivative of /(r) can be written as 

oo 

■'^ ^ T{n-a)J df^ ^ ' 

r 

where n is the smallest integer above a when it is fractional. This paper is concerned with 
the fractional derivative across the boundary and the phase (—1)""^, was chosen to make 
continuous across a — 1. Applying the derivative definition to 1/r for a < 1 (n = 1) we 

find 

D"r-^ = _r-(i+«)r(l + a), a<l (A4) 

For a — 1, this gives the usual first derivative of 1/r. For a > 1 (n = 2) the derivative is 
the same. One should be careful not to interpret the derivative for a = 2 as the second 
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derivative. The second derivative would follow from a double application of D°'. For this 
function the single derivative at a = 2 is not the same as the double application of the 
derivative operator. 

APPENDIX B: FLUID SPHERE FORMALISM 

Consider the general static spherical metric over the interior fluid with i/{r) and A(r) 

ds^ = -e'^dt' + e^dr' + r'dQ^ (Bl) 

With Einstein's field equations as Gij — SirTij, the energy-momentum components are 

SnTl = -e'\X'/r - l/r^) - l/r^ 

SttT; = e-^{iy'/r + l/r^) - l/r^ 

Stt^ = SttT; = (e-V2)K + ii^'/2 + l/r)(^' - A')] 

For the fluid interior, the energy-momentum, with four- velocity t/* = (e'^/^, 0, 0, 0), is 

T'^ = {p + P)UW + Pg'^. 

In the comoving frame the fluid stress energy is 

It is common to use the function m(r) in g^r with 

= [1 - 2m(r)/r]"\ 

so that 

A' = 2{m'/r - m/r^){l - 2m/r)-\ 
From the first field equation we have 

m' — A-KT^p 

The second field equation provides a relation between the fluid pressure, P, and 

v'/2 = {A-ktP + m/r^)(l - 2m/r)-^ 
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